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ABSTRACT 


The inireetigation of the vibrations of a tusribine blade 
has been a fasolnatlng and ehallsanglng problem to the design 
engineers and the researeh Kfosdcers and they have 
i) idealired it as a beaiB| a plate and a shell, 
ii) incorporated one after another the additional effects of 
asymmetry of the eros8«*seetion, taper, pre*^twlst, stagger 
camber, paeketting, rotary inertia, shear deformation 
eto« and 

ill) utilized differ<iHftt approximate/nimexical techniques of 
solution like Hayl«dgh^Ritz, Rltz-Qalezkin, Holzer, 

A 

Collocation, FDM, FBM etc# 

Here, the free vibration of a stationary blade has been 
dealt with modelling it as a beam and as a plate# The finite 
element technique using Gal exkin* s method has been employed to 
reduce the system of partial differ«itial equations to a matrir 
6igsn#*value problem# The ^get-lvalues and the eigen«>vectori 
were computed using a NAG subroutine FOSAEF srttich makes use of 
the Householder’s decomposition and the Qt algorithm# As the 
first eiget-value, thus computed, esdtibited oscillating trend 
imther than converging to its lower bound with the increase in 
the number of elemetts, in lieu of the eigen-value problem, an 
inverse elget-value problem was solved* 



To dotoiralno whothor a blado of asysBaotiio orosa*aeettori 
can be treated at a beam^ one neede to compute tbe chord Iwgrth 
and the maximum thickneea of the croaa^teetion* This has be«f) 
done invoking Zidcovsky*s transformation* While treating the 
blade as a beami the asynmietxy of the cxoss«*8ection and the 
pre-»t«dst accompanied by the fibre-»bendlng during torsion were 
accounted for* The freguencles were computed for three differ- 
ent data of the blade* The effect of pretwist on the coupled 
bending-bending modes and the torsional modes for a beam of a 
rectangular cross-section was studied* When treated as a plate» 
the blade was assumed to be so thin that the effect of the 
rotary inertia and the shear deformation could be neglected* 

In computation of the element stiffness matrix and the element 
mass matrix# integration was carried out by a Qauss quadrature 
formula* The frequencies were computed for four different 
aspect ratios of the blade* All the results obtained were 
compared with those obtained by the respective authors* 
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INTRODUCTION 

The turboiaachines are vdcJely used in the steam tur- 
bines,, gas turbines, compressor and gas turbines of the air- 
craft propulsion systems, marine drives etc* The steam turb- 
ines comprise of three stages called as high pressure (HP), 
intermediate pressure (IP) and low pressure (LP)t each stage 
consists of the moving blades fixed to the rotor and the stat- 
ionary blades fixed to the casing. The turbine blades are the 
soul of the whole turbomachine* The minimum satisfactory height 
of impulse blading is from 0,1 to 2*5 cm in small turbines and 
from 4*0 to 5.0 cm in lairge turbines 0^5* values of exit 

angle for the turbine blades in vogue vary from ISO® to 165® 
at HP and IP stages and from 140® to 150® at the LP stage* 
sometimes reaching 130® to 140® in the final stage of large 
turbines where maximum flow is needed* The values of the 
chord length of the bladings of the impulse turbine range from 
4.0 to 5.0 cm at the HP end to 15 cm or more at the LP end* 

The HP turbine blades are short, and tapered, have nearly equal 
entrance and exit angles thus requiring mnall angle of the 
pre-twist say of the order of 2®, have stagger angle between 
25® to 35®, and are binded at a definite pitch* The LP turbine 
blades are long wiiii considerable pre-twist* In general, a 
turbomachine blade is a tapered, pre-twisted, staggered and 
cambered cantilever of asyiomaetric alrofoil cross-section* 
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B«Gau$« df the c«aatsei<jl and the flexural eaatre ef 

the aesofoil aeetien <Jb not coincide* 

The working fluid, throttled through a noixlti flows 
over a blade and in due course exerts a tangential force on 
It, thereby causing it to rotate about the turbine axis and to 
undergo bending and torsional deformation* The process of 
power generation in the turbine leads to the excitation of the 
blades because of the flow pattern of the working fluid* 

Cyclic stresses are set in the blades due to forced vibrations* 
They are detrimental to the blades as they are likely to cause 
fatigue failure* Moreover, if the forcing freguency is close 
to the natural fregusncy of the blade# it leads to the pheno** 
menon of resonance thus further enhancing the amplitude of 
cyclic stresses and hence increasing the probability of fatigue 
failure of the blade* Therefore, it is Imperative to know the 
natural frequencies of the turfcdne blade* 

The vibration analysis of the turbine blades has been 
evolved from using the simple beam theory to that of the 
sophisticated finite element techniques* When idealised as a 
simple cantilever beam of symmetric cross-section (rectangular), 
the analysis of vibrations of the blade is quite simple* The 
analysis of the blades with aerofoil section has be«i arrived 
at, incorporating successively the effects of taper, asymmetry# 
pre-twist, fibre-bending, Coriolis force# centrifugal force, 
stagger angle, rotary inertia and shear deformation* Additional 
effects due to factors like residual stresses# relief of 
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«tr«$ges of <sarafo«r«d ggotion during bimdlng^ th« ponotration 
of tho strogg pattgxn Into tho root, the angle at the root and 
tip to the axig of the blade» non-qioincid<»it flejeoral and tor* 
gional eentreg ete# have been gtudied by vasioo# regearchere* 
Modelling the blade as a oanti lever beam with agvmmetric 
aerofoil section, pre-twist and taper can give good results 
when the length of the blade is much greater than the other 
dimengiona# However, a short fiat blade has to be idealized 
as a thin plate and a short twisted blade a thin shell* The 
effecte of rotary inertia and shear deformation have also bees 
inoorpoj^ted and studied by many researchers and scientists 
[tit, BSD, m* 

When idealized as a uniform straight beam of as-ymmetric 
aerofoil section, the blade exhibits the coupled vibrations of 
flap-wise bending, chord-wise bending and torsion* This is 
owing to the asymmetry of the section by virtue of which the 
centtoid and the flexural o^tre are non-coincident tidiereas 
wh«a considered as a beam of rectangular «rogg-*sectioii (the 
centroid and the flexural c^tre bidng coincident) coupled 
bending-bending vibrations occur which are attributed to the 
pre-twist* Because Of pre-twist the product of inertia comes 
into picture which couples the flap-wise bending and chord- 
wise bending! but the torsion remains uncoupled# The only 
factor to couple torsional modes to tho bending modes is 
distinct location of centroid and flexural centre* If the 
blade is assumed to behave as a flat rectangular cantilever 
plate, no coupling occurs between bending and torsional modes, 
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a$ thd displacements along the axis and the chord are direotXy 
proportional to the slopes along the axis and the chord res^ 
pectlvely^ Neirertheless, in the case of the plate also^ 
coupling occurs only if it is pro»twisted in which the two 
dimensional element has to be treated as a shell « In all the 
eases discussed above, it is difficult to obtain the closed 
form solution* Mence one has to resort to the approximate/ 
numerical methods like Hayleigh«-Ritr, Holrer, Collocation, 
Galeadcin, FDM, FEM etc. 

BRIEF REVIEW OF THE AVAILABLE LITERATURE t 

Houbolt and Brooks HQ developed the differential 
equations of motion for the lateral and torsional deformations 
of twisted rotating beams vdthout the loss of generality as 
regards the non*»colneidence of the torsional, flexui»l and 
centroidal axes for the application to helicopter rotor and 
propeller blades* Carnegie HQ derived the governing equations 
of coupled bending-bending*-1xjrslon motion of thin pre-*twisted 
cantilever blade using variational method* He incorporated 
the additional effect due to pro** twist while studying the 
vibrations of a straight, pre.»twisted uniform blade of rect-. 
angular cross-sec^on as a special case* Paranjpe HQ studied 
the self-*su8tained oscillations of the turbomachine blading 
lidiich were likely to occur particularly when light sections 
were subjected to heavy aerodynamic loading* Carnegie H4ll 
derived the general equations of motion of the rotating canti- 
lever blade applying the variational method* Kulkami et al HlQli 
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appiiftd both the fielener method and the total potential 
energy approaches to deduce the equations of the coupled 
bending-obending vibration* of the pre-twlated uniform rect- 
angular cantilever blade accoimtlng for the «hear-deformation 
and the rotary inertia effect* They found out the natural 
frequencie* by Ritr process and compared these with the 
experim^tal results for a specific data# 

Duggan and Slyper applied a modified c^tral 
finite dlffermoce method to the equation# of motion to compute 
the natural torsional frequ^cies of vibration of the pre- 
twisted cantilever beams of non-uniform symetric cross-section 
considering the effect of fibre-bending* Belgawakar et al QOll 
obtained the natural frequ^cles of the torsional vibration 
of a cantilever beam of rectangular cross-section vdth uni- 
form taper by collocation method and compared these vdth those 
determined by Holrer's melted LC]# Caamegte and Bao [I7j app- 
lied Qaledcln’s procedure to determine the natural frequencies 
of an untvdsted cantilever blade vdth an asymmetric aerofoil 
cross-section executing coupled b^iding-bending-torslon 
vibrations# neglecting the effect of the fibre-bending* Rao 
Q8D computed the coaled b^ding-bendlng frequ«icies of a 
pre-tvdsted cantilever blade of the uniform rectangular cross- 
section* He determined the coupled bendlng^bending- 
torsion frequenosies of vibrations of a laiiform untvdsted 
cantilever beam of an asyrnnetric aerofoil section using the 
collocation method and ignoring the effect of the fibre- 
bending* He D® computed the natural frequencies of the 
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tarsiofiai vibrations of a pre-ftwlsted cantilover boani with 
the unifoim rectangular cross-section* 

Camegie an«3 Rao flf] found out the natural frequencies 
of the torsional vibrations of the tapered* pre* twisted 
cantilever beams of the rectangular cross-section* Qupta 
and Rao S#S* [213Q used finite element method (FEM) to study 
the vibration of a single* tapered, pre-tvdsted blade of an 
aerofoil cross-section* Bajaj and Dhoopar Ql4] studied the 
vibrations of the packetted, tapered blade of the rectangular 
cross-section with the negligible pre-twist applying FEM. 

Kara and Subrahmanyam LlO solved the problem of calculating 
the natural frequencies and the mode shapes of the rotating 
blades by a finite difference method (FDM)* Saurabh Kumar [[iC! 
tackled the*probleoa of the coupled bending-bending- torsional 
vibrations of a single cantilever blade with the effects of 
the pre»twist, rotation* hub-radius, stagger angle and the 
asymmetry of the cross^-section by FBM using the Galerkin* s 
technique* Rao Lt2l] compiled the results pertaining to the 
effects of the speed of rotation* asymmetry* stagger angle 
and hub radius on the natural frequencies of the vibrations 
of the pre- twisted cantilever bladings used for the turbo- 
machines* 

In their state of art paper, teiasa* Klelb and Macbain 
Df] compiled the results of vibration of the twisted canti- 
lever plate obtained by a Joint government/industry/university 
research study* The nimerlcal results were obtained for 
tw«aty different geometries using nineteen different methods 
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15 using FEM, 2 using shell theory, 2 using beam theory — in 
order to have comprehensive, definitive set of results. The 
graph of the non-dimensionalised frequency versus the pre-twist 
exhibited different trends in some cases, so much so, that the 
predicted dependency of the frequency on the twist-angle for 
the twisted cantilever plates spanned the spectrum forom incre- 
asing significantly to decreasing significantly. This, indeed, 
made the problem quite challenging and fascinating. 

SCOPE OF THE PRESEMT WORK ; 

The present woadc aims at studying the free vibration of 
the turbine blade treating it as a) a beam and b) a plate by 
finite element method (FBvl). It is observed that the effect of 
rotation of the blade is negligible, hence we study the vibra- 
tions of non— rotating blades. To limit the study to a manage- 
able size, the effects of stagger, camber, taper, centrifugal 
force, Coriolis force, rotary inertia, shear deformation, packet- 
ting, hub-radius etc. are neglected. 

As reviewed above, vibrations of the blade treating it 
as a beam have been analysed by various authors using classical 
approximate methods. The objective of this work is to re-analyse 
some of these problems using FEf^i thereby establishing that FEJvl 
can be used sufficiently accurately to obtain the natural frequ- 
encies of straight or pre- twisted long blades by treating it as 
a beam. FEM has been used by some authors earlier also but not 
so extensively. It is also established that a short straight 
blade has to be analysed as a plate and here also FB/i can be 
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used with reasonable accuracy. A short twisted blade has to be 
modelled as a shell. Due to lack of time, this study could not 
be undertaken. 

Chapter 2 deals vdth the problem of vibration of a blade 
considering it to be: 

1) straight, untwisted beam of asymmetric aerofoil section in 
coupled bending-bending-torsional modes caused due to the 
asymmetry of the section. The data used in the reference 
[ 12 ] to solve the problem is used to compute numerical 
values of the natural frequencies and to compare the results 
with those obtained in the same paper. Coupled modes of a 
straight cantilever blade are studied assuming that the 
principal axes of the aerofoil section are closely located 
to the reference axes, say of the order of 2°. In other 
words, the product of inertia with respect to the reference 
axes is assumed to be negligible. The effect of fibre- 
bending is ignored in this case. 

2 ) Straight beam of rectangular cross-section — untwisted as 

well as pre- twisted. For untvdsted blade, flapwise bending, 
chordwise bending and torsion are all uncoupled. In the 
case of pre- twisted beam, flapwise and chordwise bendings 
get coupled due to pre-twist but torsion remains uncoupled 
as centroid and shear-centre of the section are the same. 
Using the data of the references [18] and [8] and incorpor- 
ating the effects of pre- twist and fibre-bending, the 
results are obtained by and compared with those given 

in the same paper. Using the data of the reference [17] 
and ignoring the effect of fibre-bending, the natural 
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frequencies are computed and the results are compared with 
those compiled in this paper. 

Chapter 3 consists of modelling the blade as a simple canti 
lever uniform thin plate of rectangular cross-section. The flap- 
vdse bending, chordwise bending and torsion are not coupled becaust 
displacements along axis and chord are proportional to the slopes 
along the axis and the chord respectively. In this case, the data 
from the reference [17] is used and the natural frequencies obta- 
ined are compared with those given in the same paper. 

In each case, as the system of partial differential equa- 
tions is coupled and its closed form solution is difficult to 
obtained mathematically, the system of equations is solved by Finit 
Element Method. This comprises of discretizing the blade into a 
finite number of elements, approximating the deflections over each 
element by suitable shape functions, computing element stiffness 
and mass matrices thereby converting the system of the partial 
differential equations to a matrix eigen-value problem, assembling 
the element stiffness and mass matrices, applying the essential bou 
dary conditions and solving the eigen— value problem. Finite elemen 
method, as described above, is employed because results obtained by 
this method converge to the exact solution with the refinement of 
the mesh of elements. Two noded beam elements and four noded plate 
elements have been employed in Chapters 2 and 3 respectively. To 
solve the eigen value problem, a NAG subroutine titled F02AEF [30] 
is made use of. This subroutine uses Householder's decomposition 
and QL algorithm to compute eigen-values and eigen vectors. The 
graphs of Chapter 2 have been plotted using GPGS graphic subroutine. 
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CHAPTER 2 

VIBRATIONS OF BLADE AS A CANTI LEVER BEAM 

2.1 BLADE AS A BEAM t 

Wh«n the length of the blade la large a» eorapared to 
its other dimensions, it can be treated as a beam. This 
chapter deals vsdth the problssn of free vibration of the stat- 
ionary turbine blade idealized as 

a. untwisted straight beam of asymmetric aerofoil cross- 
section and, 

b. untwisted as well as pre-tvdsted straight cantilever 
beams of rectangular cross-section. 

The effects of stagger angle, taper, rotary inertia and shear 
deformation'" are neglected. 

In the case of an asymmetric aerofoil-section, its 
centroid and the centre of flexure are distinct. This gives 
rise to coupled bending-bending-torsional vibration even vdien 
the pre-twist of the beam is z^ero. The turbine blading has 
asymmetric aerofoil section with one principal area moment 
of inertia much greater than the other. 

For an untwisted uniform beam of rectangular cross*- 
section, it is observed that the flapwlse bending, the chord- 
wise bending and the torsional modes of vibration are all 
uncoupled. But, in the case of a pre- twisted beam with rect- 
angular cross-section, bending-bending coupled motion occurs 
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due to pre-twiet and torsion is still uncoupled as the Qmt» 
roid and the flexural centre are coincident* 

2*2 DERIVATION t 

(insider a pre-tvdsted straight cantilev^er beam of an 
asyiaraetrlc cross-^section as shoum in the Figure 2*1* Let 
u, V and © be chordvdse displacement, flapwlse displacement 
and angular displacemeat about the axis of torsion of the beam 
respectively* The strain energy il'is beam is given by 

the expression: 

V - i j(eClyy(^) + 3 + 

Where I^„, I and I are the area moments and product of 

inertia of the cross-section of the beam about the reference 

axes* The term containing is because of the pre-twist r 

xy 

and it causes coupling of flapwise and chordwise bending* 

Let 1^2 principal area moments of inertia of 

the cross-section* The area moments and product of inertia 
depend on the angle betweim the reference axes and the pxin-<* 
cipal axes of the section and their transformation is given by 

I ** I ' 1 + I . A. (see referenced ) 

* av * ran *'b 


where 
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•^00# 


2ll 


av 



and I 


mu 


f.22 Zl xi 

2 


(2.1#b) 


C is the torsional stiffness of the cross-section and is glv^a 
by 

2 


r * 

C #■ + O^ij) 


(2,1*c) 


The term containing C 2 in equation (2*1*a} is the contribution 
of fibre-bending during torsion* 

G<j> Cg and Cjj are the functions of material propertita 
and geometry of the beam. The terms corresponding to these 
constants are due to the total torque on the cross-section 
that comprises of 

i) St* Vendnt torsional momentt This is because of the 
twisting as well as warping of the cross-section,, 
li) torque due to pre- twists The pre-twist augments the 

warping of the cross-section thereby increasing the 
restoring torque and 

ill) torque due to fibre-bendings The tendency of the 
beam to reduce warping thus trying to retain the 
planeness of the cross-section during the deformatlen 
results in the additional torque which is the consen* 
quence of the bending of each individual fibre* 
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where T and r are the co-ordinate distances between the 
X y 

centiold and the centre of torsion (vdiich it assumed to be 
coincident with the flexural centre of the cross-section)* 

The term containing r and r in the equation (2«2*a) is due 
to asymmetry of the cross-section w^ich is responsible for the 
coupling of torsion and flapwise as well as chordwise bending* 
If the cross-section is symmetric, say rectangular, the cent- 
roid and flexural centre are not distinct and the torsional 
mode is decoupled from the binding modes# is the polar 

moment of inertia about the centre of flexure and is given by 

Where I is polar mora^t of inertia about the centroid of 
on 

the cross-section# 

Applying Hetiilton* s principle, 

8 / (T - v)dt 0 , (2#3.a) 

One obtains the govetning equation of coupled bsnding-bending- 
torsion vibrations as 



f4 


[^A«'3”+|B” + C« (2#3*b) 

in which dach {* ) and dot (•) repraawit partial darivative# 
with respect to spatial co^^oixllnate z and time t respeetlirelyf 
and 


T 

U W V 

^ , A * 

El 

***-. yy. 

El 

xy 

0 

^xy 

^^xx 

0 

0 

0 



gA 

0 

9 Ary 

B diag \0 

N 

Oil 

1 

o 

0 

fA 

-9Arj^ 



9Ary 

-?Ar^ 



(2,3.c) 

It is obvious that tensor A causes bending-bending coupling 
whereas C causes bending-bending- torsion coupling* 

W 

At the fixed end, the deflections and the deflection 
derivatives are zero* Hence, the geometric boundary condi- 
tions are 

at z e= 0, 2 i *= 0 and u* = 0 (2*3*d) 

At free end, the shear forces, the bending moments and the 
tvdsting moment are all zero* Hence, the natural boundary 
conditions are 

atiwl, 

® H* ** £ 


(2*3*e) 



First of all, consider an ttntwtsted «trai§^t feeara of 
asYSHBotric aerofoil cross-section in which case the effect 
of fibre bending is negligible i*e*, for r « 0, Cg i» sraall 
and hence neglected and 

A «= diagpigg Oj (2*34f) 

The tensor g remains unaltered resulting In bending-bending- 
torsional coupled mode of vibration* 

As closed form solution of the system of equations 
(2*3*fo), (2*3»d) and (2#3*e) is difficult to obtain, to solve 
it, one has to apply one of the approximate ntanerical method# 
such as GaleiJcin, Collocation, FDM, FEM, etc* Since the 
results obtained by FEM converge to the exact solution with 
the refinanent of the mesh, FEM is utilised to solve this 
systsm of equations* 

For this case, since Cg is small, in the variational 
functional of the equation (2.3*a), the highest derivatives 
of u and v axe of the second order and that of 9 of the first 
order* Hence, the principal derivatives are (i) 'fj and u, 

(ii) ^ and v and (iii) 9. The compatibility condition states 
that the approximating functions should foe such that u* Ir# 
u, e should be eontlnuous aeiose the loterelement boutidasy. 
The highest derivative occurring in the variation is of the 
order three for u and v and of the order for one 9 (as stfWi 
from the natural boundary conditions, equation (2*3*e))* 

Hence, the completeness condition implies that the shape 
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function# for « and v shotild be cubic polyneralal# and these 

for 9 should be linear polyncMeials. 

The beam is discretired into a finite number of two 

noded elemwot#* A typical element is shoim in the Figure 2#2# 

au- 

There are five degree# of freedom at each node* Let sT*** 

1 

be local degrees of freedom at one of the node# 

^ av 

of the elffnent and U 2 * v-g# ' ®2 '^bose at another 

node* One can then express u at any point P on the element a# 


- ..e ia^t 
u = S U e 


(2.4*a) 


where 
(0®) • 


BU 

L«i ^ 


av 


V 


1 ^2 


e 


Oft Vft 0^ J 


2 32 


2 "Sx 




and 




«SB 


Pg 0 0 0 P3 0 0 0 

0 0 Pg 0 0 0 P 3 0 

0000 0000 Qg 


( 2* 4« b ) 


Here Pj,*s (i « 1 to 4) and Qj*s (j » 1 and 2) are the func- 
tions of the local co-o»dlinate T* 

For e elemwt, the relationship between X and the 
global co-ordinate z is given by 

^ f (* - Zq) 


where 



z 


e 


a * 
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is the half-length of the elernent. 



is the global co-ordinate of the mid-point and 
Zg and are the global co-ordinates of the nodes. 

(2.5) 

The expressions for P^’s and Q j * s which satisfy the 
convergence criteria are 

= 1 (2 - 3S + |3) 

Pj = I (1 - 5 - f * t®) 

P 3 = 1 (2 + 3S - i) 

( 2 . 6 ) 

P4 = f (- 1 - £ + + I®) 

Q, “ ^ (1 - £) 

Qj = 5 (1 + 0 

Substituting u = S U® as an approximate solution 

in the equation (2.3.b), one gets the residue 

^ = U 4 i”l " + I i" ^ y® - g s y® (2.7. a) 

£ is a non-zero vector albeit it tends to zero as the size 
of the element is reduced. 
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Applying Galpikin’s method, one gets 


•0+1 


f 

z 


sT E® dz 


(2*7*b) 


e 


whenee the system of partial differential equations simplifies 
to the following Eigen-value problem 


S ® - A i ® if 


.e 


( 2«7 # 0 } 


where * and the expressions for the el^nent stiffness 
matrix, the element mass matrix and the elmnent stress-resul- 
tant vector are as follows* 


K 

m 


e 


*e+1 

; [ (s'^)" (A)(S)" - (|'^)' (Bl(S)’]dz (2.7.d) 


8 S 


and 




Z jjL 
0 


F- <4 M")’ + I a' + (i^y (A a«)} 

LSStaS**** mm Jfm 


' 0+1 


e 


(2#7*e) 


{2*7.f) 


To evaluate the el®nent stiffness and mass matrices, 
one changes the variable from the global co-ordinate t to the 
local co-ordinate t given by the equation (2*5)* Then 

dx » a d*^ , 

9 1 3 jt 1 5 ^ 

Tz *" a ** "2 


(2.7.g) 



19 


and the limits of the integration ohange from and to 
-1 and t respeotivelyg thtis 

1 

/ f( 2 !)d 2 «s a X f(*^ at )di (2*7#h) 

e 

After integration, one obtains element stiffness matrix and 
element mass matrix given by the equations {2*7*i) and (2#7#j) 
respectively# These are svmrnetric positive definite matrices 
of the order 10 x 10, The element matrices, thus obtained, 
are assembled together to obtain the global stiffness matrix 
and the global mass matrix according to the standard procedure 
of PEM and each of them is singular. 

To solve this system of equations, one partitions “toe 
global matrices after applying the geometric boundary condi- 
tions 

at 2 ss 0, ** 9f g ' g **0, V| 0, ■ g - g l* as 0, 0^ as 0 

<2.7#lc) 

and hence one gets an eig«i-value probl^n as 

KU-AMU as 0 (2.7.1) 

mm ^ ^ ^ 

Solution of this equation gives the eigen-values s in 
number equal to the order of K or M and the corresponding 
eigen-vectors# 
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2*3.! SOtUTION l 

To solvo the equation {2»7*1) for the eigwi-values 
end idgen^veotortf the NAQ suhxoutlne F02AEF L30] has been 
utilleed. According to thi% for the df equatloije 

[a]|x| «diere [Al ntist be a jypwetric equaagi uiatsix 

and [Bl nrnat be a poaltlve definite aywraetilc^^iQuare aiatrix^ 
all the elgen-valtijia and eigen*veGtora are <§^puted using 
B»uii^older*s decomposition and QI algorithm# Note that the 
eigen-*values \ represaait the sjtuare of the natural fregucn- 
eles 0)4 of vibration of the bews* As the con;U.nuous system is 
approjdmated by a set of fl^^ite number of el aments 4 the Infinite 
number of degrees of freedom corresponding to j:he continuous 
system is reduced to a finite niaaber, thereby making the 
approximated model stiffer than the actual aystem. .Obviously^ 
adth the increase.. in the number of elemeniSt the number of 
degreet.,.of frj|ed<wa Increases^the overall stiffness of the 
epproximated model decreasee and consaqumitly the eigan**values 

converge to their lowj^r bounds* 

During the cojgputation of ■Wie resultSf it was expert 
lanced that as the number of el.f»»ents was increased* the first 
eigen-value oscillated about a maan valua whereas the higher 
order elgen-valuet converged. It was coniecturad that this 
could possibly be due to 

(i) ixx» leaang to the ill-conditioning of the 

•ystam„,of equations f and 

(11) tha number of itorationa in the NAG aubroutine are 
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fixed but tmknewn end the wter cannot exerclee any 
eontrol <Wk'^ it* The eig«i-val«es are the roott of 
the pol'yno»^al flvan by 

I A - ^ b1 i» 0 (2*7*!h) 

The higher root* of the polynomial equation (2*7*m.X 
can be obtained by neglecting the lower degree terms 
while for finding out the lower eigen-values^ in part- 
icular the lowest one* one has to consider all the 
terms* C3onsegu«atly, higher roots can be computed 
with fewer nisnber of iteratilns whereas the first few 
roots can be obtained only after a considerable number 
of iterations* Therefore# it it qtiite apt that the 
value of the first ^$ig«ri-valuc is not a converged 
value within the number of iteration# fixed by the 
subroutine it%slf* 

Since th^t ntJBiber of Iterations is beyond the contwl 
of the user* to mltlgste the problem first 

natural frecBi***®^* solves the^^systsra sjt sn inverse elg«i- 
value problem so^^that the lowest mode becomes the highest mode 
of the new system. Thus* in lieu of solving the systws of 
equations (2*7*1)* one solves 



where M i 


^ H 

_J 

^(n-i+1 ) 


(2«7*n) 


The ncmfoer of eigen-values ^ for the new systmii of equations 
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i« ACfuftl to the order of K or The eigwtt^-values are the 
inveree of the square of natural frequency Hqnce, ... 
the flret (lowest) eiqen-value of the actual sys ten,, becomes 

the last (highest) eigen-value > 14 ^^ of the inverse system 
(yL/^ • converged value of is the reciprocal of 

th. eoniergad value It la verified that thla indeed 

yields^, converged eig.^— values with the increase in the number 


of elem^ts* The number of frequencies and the corresponding 
eigwi— vectors ..obtained is equal to the order of the partitioned 
stiffness or mass matrix and can be found out..by subtracting 


thj| number of the 8pecified...degrees of free^,om f3»m the total 
number of degrees of freedom of the approximated model* 


2*3*2 DATA AND RESULTS » 

The results have been obtained for the blade used in 


112] so that the results can be 
for the aerdfoil section are* 


4»*. 


^11 

m 

3*49634 X 

^22 

m 

2.79291 X 10**^ 

1 

m 

0*1524 m 

Xx 

m 

1*9304 X m 

r 

y 

m 

1*1938 X 10*^ m 

*♦4 

A 

m 

-/v-5 2 

5.8967 X 10 m 


compared with it* The data 


Geometric 

properties 
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E » 213*7380 GPa 

***i 

C » 9.1485978 H/m 
^ » 7858.3702 kg/ra^ 


l^at axial 
pxo parti ea 


To tast whather the given cioas-section can be 
treated t$ a beam or not, tha Zukovsky transformation 22 was 
made use of. The Zukovsky transformation, which maps an aero- 
foil-^section in ■.•plane to a circle in z-plane as shown in tha 
Figure 2.3, is given by 

2 

^ ^ , where ^ + i h 

* ' (2.8. a) 

and 2 « X + iy 


Refer to Appendix I fpr details* 

*.* “^he' transfonnation gives the chord length b and the 

maximum thickness h as 

b « (2.8.b) 

and 

h « 3(2I^^/A)^/^ (2.8. c) 

#*4 

From egoations (2.8«b) and (2*8*e) one gets 
^ *» 5.53606 

and 

I * l^^22^^11^^^^ * 8*426461 
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As ths X(mg(|h^is Xargs ia eompaiison with the Ghe3»i 
Xtisigth and the maximw! thicfcaess,^^/or the data considered, 
the blade can be treated as a beam. 

flesaXtst For the data givsn above, the computed results 
upto fifth mode are given in the Table 2*1 and are ccmipared 
S^th those obtained in the re,teren'ce Cl2l| by (a) analytical 
method and (b) Ritz-^Saleiicln method. The physical identifi- 
catioi^ of the frequencies was done by drawing the correspon- 
ding mode shapes (they arjg not shown here). The results were 
obtained for fourteen elements of equal length* 


Table 2.1 

1 b 

Comparison of results for g » 5.53606, » 8.426461 


Method 


r * *’**"' " ’ "" ’ """"" " " * — 

DescriptieaJ Frequaicy in Hz 

of t" r - r ’ ""T r * "" ’ 

modes *First {Second {First {First {Third 

{binding {bending {bending {torsion {bending 
{in YY {in YY {in XX { {in YY 

{I {II { III ? IV {V 

I x.— n, ■„»; ,,4- ...I 


Analytical results [12] 

Results by Mtz- 
Galetkiii method [1^ 


96.9 

606.5 

841.2 

1072.9 

1699*0 

96.5 

604.9 

842.0 

1089.4 

1 694.9 

96.8 

606*3 

843* 3 

1078.5 

1697.8 


Present results 


71 


2.4 CASE 2: 


Having studied the vibrations of untv\dsted cantilever 
beam of aerofoil cross-section, we now focus our attention on 
the pre-twisted cantilever beams of rectangular cross-section. 
For a rectangular section, the centroid and the flexural 
centre are coincident. Sc, the torsional mode is uncoupled; 
but the flapwise bending and chordwdse bending are coupled 
because of the pre-tudst. If the pre-twdst is zero, both the 
bending modes also get decoupled. In this case, it is assumed 
that the effect of fibre-bending is not negligible. 

If b and h are the breadth and the height of the beam 
of rectangular cross-section, the expressions for the constants 
and in equations (2.1. a) and (2.1.c) are given as 
follows: 


(i) 


n 


1 Gbh^Cl 


122 h „ 1 

5 t> ^ 

It n=1,3,5... n 


tanh 


n-n b -1 
2 h ^ 


(2.9. a) 

where is the St. Venant torsional stiffness [l9]. 


(ii) 


=3 



(2.9.b) 


where is ihe torsional stiffness due to pre-twist 

[lOJ . 


(iii) 


Eb\^ 

2 - 144 “ 


(2.9. c) 
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wher« C 2 is the constant for the fibre«-bending [12® • 

In this case, the tensors A and g of the eqtiatioii 
(2«3»b) are givai by 


#■ 


A • and C • dlag [_?A, ?A, (2*10.a) 

where 


and 



008 0 
sin ■’••cos 0 


{ 2* tO#b) 


L ® ® 

As the fibre*bending is considered, the ffiaxiisom deri- 
vative of © in the functional is of second order* Hmtee the 
principal derivatives for © are ^ and ©♦ The convergence 
critex|,a states that © should also be expressed by Hermitlan 



polynomials* Pj.*s* Hence, in equation (2.4*a), one gets 

(u®)^ L ^ ^ 


and 


$ » 

m 


„ !!!' 

2 ^x 


'1 &x 


W ''o ^ ®0 

r T £n TiTtT .-Trtrrg' ] 

r\ C'n J 


2 ^x 


2 dX 

(2*1t#t) 


P, Pj 0 0 0 0 P3 0 0 0 0 

0 0 Pg 0 0 0 0 Pg P4 0 0 

0 0 0 0 Pg 0 0 0 0 P3 P4 


(2*t1*b) 
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The expretslofis for el«Dent stiffne## matrix and mtsa matrix 
are the same as those in the previoxis case* *ith the incor- 
poration of fibre-hwieEng, the order of the eionimt stiffness 
and mass matrices hecomes 12 x 12 and can be ewalaated by 
simplifying integrations^^are discussed in the Section 2*3# 
Applying the geometric boundary conditional, 

2>V- 

at z *5 0| tJ| » 0, * Oj » 0 ,, «* *" Oj » 0 

{2*1t*c) 

which physically mean that the deflection and the deflection 
deiivatifefes are zero at the fixed .end of the cantilever beam, 
one obtains the eigei>^value problem on the lines similar to 
the Section 2*3. 

f t , 

2*4*1 DATA AND RESULTS t 

The results have been obtained for the blades used in 
D® t CCl and Hit] so that the results can be compared with 
these cases. 


2*4* 1*a Blade 1i The effect of the fibre-bending i 


Data I b m 0*0254 m 


h * 0*001727 ra 

Geometric properties 

1 •* 0.1524 m 
r » J rad 
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E = 206.85 GPa 

G = 82.74 GPa Material properties 

p = 7858.3702 kg/m^ 

"I 

As ^ = 6 and ^ = 14.706, the length is larger than the other 
dimensions and hence the blade can be treated as a beam. 

Results: For this data, the frequencies obtained are given in 

the Table 2.2 and their modes are identified as discussed in 
the previous case. The beam was discretized into 14 elements. 


Table 2« 2 

Frequencies for ^ ^ 14.706 and Y = ^ I'ad 


Mode 

No. 


Frequency in Hz 


— j — -r 

Bending in Y-Y ! Bending in X-X , 
direction (coupled); direction (coupled). 


Torsion (uncoupled) 


;Present ,'Results of {Present 
results {ref. [l8j {results 


{Results of {Present 
{ref. [18J {results 


{Results pf 

{ref. 11 8 J 


1 62.027 62.17 

2 310.005 310.18 

3 1232.751 1232.54 

4 2223.903 2223.51 


977.620 977.62 756.421 758.15 

5930.315 5926.49 2298.804 2349.61 

39^.691 4162.77 

5702.135 6242.28 


It is observed that the frequencies of bending-bending 
modes tally with the results given in the reference [18], but 
the torsional frequencies do not agree with the results. This 
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is because in the reference [lBj, 

(i) the section has been aporoximated to the elliptic 

form and the corresoonding formula for the St. Venant 
torsional stiffness has been used and 
(ii ) the contribution of the pre-tmst to the torsional 

stiffness is neglected and also it is assumed that the 
fibre-bending effect is negligible. 

2.4.1.b Blade 2: The effect of the pre-tv\dst: 

Let us now study the effect of pre- twist on the vibra- 
tions of the beam with rectangular cross-section. For this, 
the data used in the reference [S] is considered. Using this 
data, the frequencies are found out for different values of 
pre-twist* The graph of reference [s] is shown in Figure 2. 4. a. 

Data: b = 0.0254 m 

h = 0.00635 m Geometric properties 
1 = 0.3048 m 

E = 206.85 GPa 

G = 82.74 GPa Material properties 

3 

p = 7858.7302 kg/m 

For this data, ^ = 12 and | = 4 and hence the -blade is treated 
as a beam. The frequencies for vailous values of pre-twist 
are tabulated with their modes identified and their trend 
with respect to the pre- twist in the Table 2.3. The beam is 
discretised into fifteen elements. 



TaW« 2,3 

The effect of |>re- twist on the frequencies for 
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It is observed that each f^ode has got a specific 
trend of variation with respect to the angle-of-pre- twist. 

The first frequency corresponding to the first Y-Y bending 
mode increases with pre- twist. The second mode corresponding 
to the first X-X bending mode decreases with pre-tvdst and 
the fifth rr.ode corresponding to the first torsional mode 
increases with the pre- twist. This is because the beam 
becomes stiffer for the Y-Y bending and the torsion whereas 
it becomes more flexible for the X-X bending at their first 
modes. The graphs in the Figure 2.4 depict the trends of the 
first nine freouencies with respect to the angle of pre-tv.lst. 
These trends tally with those illustrated graphically in the 
reference [8], Figure 2. 4. a. 

2.4.1.C Blade 3: The effect of the pre-twist and the 

dimensions: 

For this blade, the data from the reference [17] is 
utilized. 

Data: h = 0.025 m Thickness 

E = 200.00 GPa 

VI = 0.3 (Poisson's ratio) f^aterial properties 
P = 7800 kg/m^ 

For this data and for four different geometries, the non- 
dimensionalized frequencies given by 
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where 


^Ph/D 


12(1 - |i^) 


( 2 . 12 ) 


are computea for the different angles of pre-tmst and are 

tabulated in the Tables 2.4, 2.5, 2.6 and 2.7. The beams are 
discretised into seven elements. 

Table 2.4 
F - h = ^ 


Pre- 

twist 

Y 

rad 


Y-Y 

ben- 

ding 


Tor- 

sion 


Non-dirtiensionali zed frequency 

III I IV ; V I VI I VII I VIII 

\ 1 1 j j 

X-X ; Tor- ; Y-Y ! Tor- I Tor- } Y-Y 

ben- { sion } ben- I sion } sion J ben- 
ding ; ; ding j | ; ding 


0 3.354 5.940 16.771 18.739 21.044 34.040 49.225 59.311 

11/12 3.36 3 6.154 1 6.087 18.615 21.98 9 31.5 34 45.221 58 . 487 

71/6 3.37 1 6.410 1 4.743 1 9.38 9 24.006 32.8 46 47.10 3 56.908 

11 /4 3.379 6.677 13.432 20.196 26.440 34.213 49.062 54.722 

n/3 3.397 6.954 12.257 21.034 29.086 35.633 51.099 52.315 


n/2 3.448 7.539 10.363 22.804 34.383 38.632 55.399 48.041 
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Tabl® 2*5 

1 - 1 , 20 



0 3*354 5.364 19,693 21.030 34,889 52.348 59.067 67,082 


Tr/12 3.466 9,453123.592 20.560 48.436 69.459 58,375 70*140 
V6 3.464 11,949 36,142 19.079 61,226 87.800 57.282 77.227 


Tr/4 3.395 14.167 42,850 17.211 72.590 104.097 56,123 86.339 


^/3 3.417 16,229 49.090 15,462 83.160 119,255 54*638 96,053 
V2 3*474 20.105 50,811 12,567 103.016 147.730 50.695 109.752 


Table 2.6 



Pre- 

twigi'fe 

r 

rad 



0 3.354 16,771 17,780 21,030 55.101 59.067 97*622 105,150 

TT/12 3,363 16,057 17*808 21.989 5 3,365 58,487 91.250 106,081 

3.371 14.743 17,898 24,006 54.138 56,908 91,712 108,357 

•^4 3,379 13.430 17,996 26.440 54.432 54.722 92.210 110,079 

■73 3.397 12.257 18,100 29.086 54,748 52,315 92,746 109*734 

"/2 3.449 10.363 18.333 34.383 55,453 48,041 93,939 1 06,366 
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Table 2.7 
F ^ 20 


Pre- 

twiat 

r 

rad 


Non-dim ensionalised frequency 


TT 
I 
t 

1 Y*Y 
} ben- 
} ding 

JL 


"r" ' "" *T' ■ --r ' " 

{ III { IV ? V 

■j \ ^ \ 

\ Tor- \ Tor- 1 Tor- 

{ sion { sion \ aion 

I f » 

t t I 

J I 


nr 

f 

■4 f- 

! Y-Y ! 
i ben— I 
} ding } 

J. L. 


II 


VI 


VII 


VIII 


d 

I Tor- 
} aion 


Tor- 

sion 


Y-Y 

bin- 

ding 


0 

1.903 

3.354 

5.908 

10,466 

15.695 

4f— 

9 

o 

CM 

21.030 

59.067 

Yl2 

7.269 

3.456 

21 .986 

37.246 

53.412 

70.821 

20,560 

58.375 

tr/e 

10.308 

3.463 

31.179 

52.819 

75.744 

100.433 

19.079 

57.282 

V4 

12.813 

3,398 

38,756 

65*654 

94.151 

124.840 

17.211 

56.124 

V3 

15,062 

3.420 

45,560 77.180 

110,680 

145,756 

15.462 

54.638 

Tr /2 

19.175 

3.467 

57.999 

98.253 

140.898 

186.824 

12.566 

50.595 


The trends of the frequencies versus the pre- twist for 
the four different aspects ratios are plotted in the graphs ■— 
Figures 2.5, 2.6, 2*7 and 2.8 — corresponding to the Tablet 
2*4, 2.5, 2.6 and 2.7 respectively* As surmised, some modes 
^ther exhibit the trends contradictory to those in the refer* 
smce [^tf]or do not exhibit a definite trend, for the short blade 
is treated as a beam here. 











x,y,z, = Reference axesi1,2= principal axes 
0O|= Flexural axis iPQ =Centroidal axis 
lf= Total angle of pretwist. 

FIG.2.1.a PRE-TWISTED CANTILEVER BEAM 


cfi centre of flexure; cm h centre of mass 
FIG2.1.b CROSS-SECTION OF THE BEAM 
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N B- The encircled numbers represent 
the element numbers. 

FIG. 2.2a DISCRETIZATION OF THE CANTILEVER 
BEAM INTO A FINITE NUMBER OF TWO-NODED 
ELEMENTS. 



eth element 

Z -Global co-ordinate, I -Local co-ordinate 
a -Half the length of the element 
Numbers inscribed in the triangle represent 
local node numbers. 

FIG. 2.2. b A TYPICAL TWO-NODED ELEMENT 
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FREQUENCY VERSUS PRE-TWIST 
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NON-DIMEINSIONALLIZED FREQUENCY 
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FIG. 2.7a GRAPH OF 

FREQUEMJY VERSUS PRE-TWIST FOR 
L/B « 3 AND B/H = 5 
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CHAPTER 3 


VIBRATIONS OF BLADE AS A CANTILEVER PLATE 
3#1 BIADS AS A PLATE t 

When the length of the blade is short and is compa- 
rable to its cross-sectional dimensions, it t^ids to behave 
like a plate rather than a beam* If the blades are pre- 
twisted, they have geometric curvatures and torsion and hence 
behave like shells* When their pre-tvdst is zezp, their 
geometric curvatures and torsion are zero and can be treated 
as a plate* A plate is defined as an object bounded by two 
flat surfaces# If the distance between these surfaces is 
very small as compared to l^he other dimensions of plate, it 
is called a thin plate, else it is said to be a thick plate# 
For the plates, vdnere the thickness is very small as compared 
to either the overall dimension or to the wavelmgth corres- 
ponding to the highest frequency of interest# the effects 
due to shear deformation and rotary inertia can be neglected 
vdthout affecting the accuracy of the solution DCI * 

This chapter deals with the study of free vibrations 
of a short blade treated as a uniform thin cantilever plate 
of rectangular cross-section, using FEW and neglecting rotary 
inertia and shear deformation* 

3*2 DSRIVATION t 

For the Ijiin plate shovm in Figure 3*1, let w(x, y) 
be the displacement of a point on the middle surface along 
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2:»axi«* If wa assume that the no»nal to the middle surface 
continues to be normal to it even after the deformation* 
then the displacements u and v of any arbitrary point 
(x, Y# z) are given by 


u *c z\7w i.e* u s= z 


V = z 


dw 

dw 

ay 


{ 3« 1 # a ) 


The most significant strain and stress components are 


2 

£ *= zVV w i*e* € v-v “ ^ 

2 

^ , '3 w 

^YY 9y2 

•:i2 

^ xy 9x3 y 


(3.1. b) 


and 


H I i I — ^ 2 


5C^2v w4 + (i •//) z vv*D 

•A 2 ^2 

W . n 9 W\ 


1 

l.e. s'™ - 71^ 

,2 J> 


XX 


'.y - 

^ . Bz ^ J!®a- 

(Tjjy ^ 3x 3y 


(3.1.C) 


The stress- resultant over the thickness are defined 


by the relations 
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h/2 

-h/2 


h/2 

K * S z dz 
^ -h/2 


T7 


h/2 

Ky *“ ^ iS’vv ^ 

-h/2 

h/2 

% ** ^ <fy 

^ -h/2 


Q ^ » 
Y 


h/2 

r tf*™ <Jz 
-h/2 


(3.2.a) 


The stress-jfesultant represents the bending moment per 
unit length on x-face about y-«xis. Similarly also repre- 
sents the bending memest per unit length but on y-face about 
x-tfel?#* The quantity M denotes the twisting moment per 
unit length about x-axis on x-face (or about y-axls on 
y-faee|* The shear forces per unit length along z-axls on 
X and y-faces are r<^resented by and Q^* The relation 
betwesBh the moment tensor 


¥ 

•m 




¥ M 

X xy 


“y 


and the shear force vector 


a 


r -I 

Q 


)BS 


1 
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is givsaa by 


The expressions for ^ and £ in terms of w arei 



lain 

*= D[/iv^ ^ 4 ~ /^) w) 

(3.2*b) 

and 

s 

*= D V(v^ 

w) 

(3.2.c) 

where 

n 

Eh^ 

TT" 

(3.2.d) 



12(1 



D is called the flexural rigidity of the plate* 

The strain energy [^16J of the plate is given by 


V « I J [(1 w).(vv w) +;;(V^ w)^JdA 

(A. 

tr r 2 ^ 

*= r **■ 

3x e>y ox dy 

The expression for the kinetic energy [2€3 of the plate is 
given as 

T *= 2 (3*4) 

/Applying Hamilton* s principle, 

^2 

S f (T - V) dt * 0 (3*5*a) 

t. 


one obtains the governing equation of motion of the plate as 
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D w + 




{ 3*5*b) 


At thft fixed end of the cantilever plate, the displacement 
and the displaceirient gradient are zero* Hence, the geome- 
tric boundary conditions at x * 0 are 

1 ) w w 0 and « 0 

o X 

called the essential boundary conditions and (3*5*c) 

li) IH « 0 and « 0 

li ; 'dy 21 X a y 


called the implied boundary conditions. 

At the free boundary, the bending moment, the tvidsting uMsment 
and the shear force, are all zero. That is, n = 0 and 
n « 0, where ft is the unit outward normal to the boundary* 
Hence, the natural boundary conditions for a rectangular 
plate shown in Figure 3.1 are 


i) 


At X n 1, 


oy 


^ 0 

BX3 y 


0 




— 

3 x Sy 


0 


b 


il ) At y » + 5 ♦ 


"5 A*- 2 

9y 
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a^w 

I Jjgf 

3 y 


0 


d\ 

" "" ’' VI ' " , '"’. .. ""’ + —— S 

ax^0y ^y-^ 


0 


(3*5,d) 


At the governing equation is of the fourth order, 
two boundary conditions are needed at each boundary, whereas, 
in equation (3*5*d) three boundary conditions are obtained at 
each free boundary* In order to make the boundary conditions 
compatible, the boundary conditions for the shear force and 
the twisting moment are clubbed together hence, the 

equivalent natural boundary conditions at the free edges of 
the rectangular plate are 


i) At X « 1, 


>2 >2 
L^+»2-S _ 0 

3x^ 3y 


^3 

— -w + ( 2 "• /I ) 
c>x 


3 \ 
dx ay^ 


0 


ii ) At y » ± I » 




(2 




IB 


0 


C 3«5# ©) 
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3*3 FINITE ELEiMENT FORMULATIOM : 

The closed foim solution for a rectangular plate can 
be obtained explicitly if and only if the edges of the plate 
are simply supported. However, in the case of cantilever 
plate, closed form solution is difficult to obtain. There- 
fore, to study the vibrations of a cantilever plate, one has 
to make use of some numerical technique for the given plate(s) 
and here, the FE^VI has been opted for. 


Consider the expression for the strain energy in 
equation (3.3). The highest derivatives of w in the expres- 

^ r\ 


2 2 

8 w 8 w 


sion of the strain energy are — ^ and . There 

3w ^ , 3 w 


fore, the principal derivatives are w, and The 

compatibility condition states that these principal deriva- 
tives should be continuous across the inter-element boundaries. 
The highest derivative of w in boundary condition is of third 
order. Completeness condition states that the shape func- 
tions should be such that, as the element size reouces, the 
function and its derivatives upto the highest order accuracy 
in the variation of the strain energy are of finite values. 

This leads to the conclusion that the shape functions should 
be the product of at least cubic polynomials in § and t). 

It is noted that a mesh of rectangular elements is used 
to suit the convenience of the geometry and because of the ease 
of finding the compatible shape functions. Consider a four 
noded rectangular plate element as shown in Figure 3.2. If 
(x V ) is the global mid-point of the element, p and q 


/ ^^935 
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are half the length and breadth of the eluent re 8 pectively» 
then a point within the element represented by the global 
co-ordinates (x, y) can be expressed in terras of the local 
co-ordinates as 



where -1 ^ < 1 

and -W ^ < 1 

The displacement of any point within the element in 
terms of the nodal degrees of freedom and the shape functions 
is 


w = S . W® 


(3.7) 


where the nodal displacement vector, is given by 




9'^^^ A 2 

3x ay 3xa y ^^2 c)x ay " ay 


aw. . -aV 

yf , , . ^ ^ „ ^ 

3 ■ax ay 3xsy 4 ax ay 'axciy 


C 3*s • a ) 


Thus, there are sixteen degrees of freedom per element* The 
shape function vector is, 

S b: 0^ 02 03 04 05 0^ 07 08 09 0|o 0^^ 0-(2 ^13 ^t4 ^15 ^16 

( 3 . 8 .b) 

The expressions for 0^’s which satisfy the convergence crit- 


eria are: 
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For 1 = 1,2, 3, 4, 

i 

^4i-3 “ 16^2 + 3X^k (2 + ^ - 7 ^^^) 

^4i-2 *= -1 37 ^ 7 - 7 ^ 7 ^) 

%..1 *= ^ 1^2 + SC^f - C^-C^) (- 7 ^ - 7 + 4 7 ^) 

^4i “ - ■§ + liT^ + f- 7i • 7 + 7i7^ +7^) 

(3.9) 

Substituting w = s . W® in equation (3*5), 

one gets the residue 8 ® 


e® = D( s ) . w® - § h s . w® 

wkicl^ cav) loe ve. Written as 

f e 


, “K 4 

3 X ^ '9 y' 


4 

^7 6 

— !5 • w 4 2( 1 - 


- §ha>^ s . w® 


3 ^s 

A)( : ' ”‘^" ' f7 2 ^* 

3 X o> y 


( 3« 10*a) 


Applying Galeikin’s method, 


Vi Vt 

i / S fi dx dy s= 0 

Y "” 

e 


Xe+I r r 3‘*S 

^ ^ l'^ -L^ " 


3 V 

^ + 2M 

5y' 




5 +^- 5 + 2(1 


- Shcj^ S Sdx dy)'™® = 0 


(3.10*b) 


Integrating by parts, the partial differential equation 
reduces to 
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K® W® = R® 


(3.10.C) 


® SS ^ 2S are element stiffness matrix;, element mass 

matrix and element stress-resultant vector and are given by 


K® « D r " r’ ’■! (1 s).(vv s) + 
^e Ye 


*9+1 ^e+1. 


H sXv^S) \ dx dy 


„ frr fs 


a^s "S^s T 

•" -'"'fTayaSTyr'^’' 


(3.10,d) 


M 


e 


9 h JJ S S dx dy 


( 3* 1 0» e} 


R 


^ ^ 04*1 ^2 ^ 

* ^ I S( -— " g + ■2)( s • W ) — ^( — -g +j^ — g)(S. W®) 


'e 


3 X 3 X 5 y' 




^ / e ~I 

< 1 < 2 - 2 ")J, dy 


X 


e+l 


+ D T s(^ + £ gXs. \f) - 's=(~-K +/W ^)(S, w®) 

tfy 3 y Sx ^ 37 “^ ^x 


' avN . 

y 
y. 


^2 


eT'®+1 

- -^)a-| 5T7 'S-E®)]^ dx 


(3.1^.f) 


The element stiffness matrix and elauent mass-matrix 
are computed Invoking Gauss-quadrature formula. 


% 
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; i fix, y) dx dy « 2 : 2 : W.w. fix,, y . ) (3*1t) 

1 -1 is1 j=:f ^ J i J 


in which Wj^*s are w^ght functions corresponding to the Gauss 
points (X|^t Vj), ngx and ngy bd.ng nrabers of Gauss points 
along x-direction and y-direction respectively# In order to 
have the limits of integration -1 to 1 as indicated in equa- 
tion (3*11), one transforms the variables in (3»10.a) and 
(3*10#b) into natural co-ordinates using the transformation 
given by equation (3*6). From this transformation, one obtains 

dx dy = pq 

(3.12) 

I 

d 

'^X'dY Pg dtdfji 

Therefore, the element stiffness matrix and element mass matrix 
are obtained as 






JL ^ 

P'1 3^2 5g2 


q 




2(l=m 








■I 


did^ (3.12#a) 


and 
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1 

^ h pq ^ 
-1 


J dld^ 


(3.t2*b) 


As the limits of integrations in the equations 
{3*12*a) and (3*12*b) are in accoitlance with those in the 
equation (3.tl), the integrations can be evaluated by Gauss* 
quadrature formula* The elsnent stiffness matrix and mass 
matrix, are svmmetric positive definite and are of the order 
16 X 16* The elem^t matrices of all the elements are assem- 
bled according to the standard procedure leading to the global 
stiffness and mass matrices which are partitioned after speci- 
fying the essential and the implied geometric boundary condi- 
tions from equation (3*5 *3) so as to obtain an eigen-value 
probl^n 

r K - ^ m 1 W « 0 (3.13*a) 

L » « J ^ 


To solve this system, i*e. to obtain eigen-values and 
eigen-vectors, the NAG subroutine titled F02AEF is used# Its 
details have been discussed in Chapter 2. As in the case of 
the beam, here also the inverse eigen-value problan is solved 
considering the system of equations to be 

(3*13*b) 


« 0 where p, 


3*4 DATA AND HESUtTS t 

The results have been obtainad for the blade consi- 
dered in reference The data is as follows? 
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h = 0.025 m 

R = 0.3 

E = 200.0 GPa 

P = 7800.0 kg/m^ 

For this cata, the non-dimensionalised freauencies given by 

o' ph X ^ 2 

“o = (3.14) 

are computed for four different geometries. The first eight 
non-dimensionalised freouencies corresponding to each of the 
four geometries are given in the Table 3.1 below. 

Table 3.1 

Frequencies of the plate 
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of nodes. The results are fairly close to those given in the 
reference [17] but do not exactly tally with it because the 
convergence of the frequencies with the increase in the 
number of elements is comparatively slow in this case and 
hence it requires more number of elements to acquire the 
desired accuracy. The increase in the nuimber of elements is 
not possible owing to the limitations of the computer memory 
space required to store five square matrices and a few arrays. 
Needless to say, more accurate results cannot be obtained. 


/ 



ABCD= Neutral plane 
Q=(x,y^o) P5(x.y,z) 
w= displacement of Q 
u.v = displacements of P 
along x and Ydirections. 


FIG. 3.1 FLAT PLATE 
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r'x (nx+1) 


nx= number of elements in the x direction 
ny=. number of elements in the y direction 
N.B. The encircled numbers represent the element numbers 

FIG-3.2.a DISCRETIZATION OF A CANTILEVER PLATE INTO 
A FINITE NUMBER OF FOUR-NODED ELEMENTS 

i P 


A (-1^1) 



p^q- half the element length and breadth^ respectively. 
N.B. The numbers inscribed in the triangle respresent local 
node numbers 

FIG.3.2.b A TYPICAL FOUR-NODED PLATE ELEMENT 
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CHAPTER 4 
CONCLUSION 


4.1 BLADE AS A BEAM t 

Table 2*1 shovfs the natural frequencies corresponding 
to the coupled bending-bending-torsion modes of a long *= 
5*53606), untwisted, cantilever beam of asymmetric aerofoil 
cross-section* The results are in agreement with analytical 

results of the reference [12l* 

Table 2.2 shows the natural frequencies corresponding 
to the coupled bending-baadlng modes and uncoupled torsional 
modes of a long {| * 12), pre-tvdsted and cantilever beam of 
the rectangular cross-section* The results are in agreement 
with those in the reference [1^ • The frequencies correspon- 
ding to flapwise bending and torsion increase w4»tb the 
increase in the angle-of-pre-tvdst while those corresponding 
to chordwise bending decrease* This is consistent with our 
expectations, as the area moment of inertia about X-X and 
the product of inertia of the section increase, the stiffness 
due to pre-twist and the effect of fibre-bending increase 
while the area moment of inertia about Y-Y decreases* ^ 

Table 2.3 also shows the results for a long beam « 

12) and there is a definite trend as far as the variation of 

the frequency with respect to the pre-twist is 

on the other hand, the results of Table 2.4 are for 

a ..oxt baa* (| - 1 « 3). Har. so*, toxalonal ox flap.1.. 



bendina frequencies decrease, chordwise bending frequencies 
increase while some frequencies do not exhibit any definite 
trend. This leads to the conclusion that 
(i) Long blaoes can be analysed sufficiently accurately 
using a beam theory but not the short blades. 

(ii ) For long blades there seems to be a definite trend 
for the variation of frequency with respect to pre- 
twist. 

4.2 BLADE AS A PLATF i 

Table 3.1 shows the natural frequencies of a short 
blade for various length to breadth and breadth to thicknes 
ratios. These frequencies are fairly close to those given 
in the reference [l7] . It is observed that the frequencies 

jL 

for the same ^ ratio are the same irrespective of the 
ratio. One can, therefore, draw the inference that 
(i) Short flat blades can be analysed satisfactorily 
using the plate theory, and 

(ii) The frequencies of short blades are independent of 
^ ratio as long as it is sufficiently small. 



63 


REFERENCES 


1* HOUBOLT JOHN C# and BROOKS GEORGE W. , 

"Differential equations of motion for combined flapwise 
bending, chordwise bending and torsion of twisted non- 
uniform rotor blades", 

NflkCA Technical Notes, 1956# 


I2)i CARNEGIE W# 

"Vibrations *of pre-tvdsted cantilever bladingt An addi- 
tional effect due to torsion". 

Proceedings of the Institution of Mechanical Engineers, 
pp* 315-322, Vol. 176, No. 13, 1962. 

3. PARANJPE mkmD A., , „ 

"Self— excited vibrations of turbomachine blading , 
Journal of Applied Mathematics and Physics IZAMP), 
pp. 598-608, Vol. 14, Fasc 6, 1963. 


6. 


CARNEGIE W., ^ ^ ^ 

"The application of the variational method to ^ 

the equations of motion of vibrating cantilever blading 
xo iirQ on ^ I 

Bulletin of Mechanical Engineering Education, 

Pergathon Press, pp. 29-38, Vol. 6, 1967. 

DUGGAN A. P* and SLYPER H.A. , 

"Torsional vibrations of pretwisted cantilever beams , 
International Journal of Mechanical Science, 

Pergamon Press, pp. 871—883, Vol. 11, 1969# 

BELGAUMKAR B.M., CARNEGIE W. and RAO J.S., ^ 

"Torsional vibration of a cantilever beam of rectangular 
cross-section with uniform taper". 

Bulletin of Mechanical ^gineering Education, 

Pergianon Press, pp* 61-67, Vol* 9, 1970* 

"Solution^of^the^quations of motion 
bending torsion vibrations of turbine blades by the 
method of Rit«-Gale3*in", c^4 

International ’ 

Pergamon Press, pp* 875-882, Vol* 12, 1970. 


8 , 


vibration of pra-twiat«J beam, of ractangular 

jJSmal'of the’ , 971“**^® ’ 

pp, 62-64, Vol* 23, No* 1, February 1971. 



64 


9. 

^Couplad bmding fotndirig torsion vibrations of canti- 
lavtr beams”. 

Journal of the Aeronautical Society of India, 
pp* 265-268, Voi* 24, No. 1, February 1972. 

10. mo J.S. , 

“Torsional vibrations of pre-tvd.sted cantilever beams". 
Journal of the Institution of Engineers (India), 
pp, 211-212, Vol, 52, No, 9, Pt. CI5, 1972, 

11. CARNEGIE w* and HAO J.S., 

"Torsional vibration of pre-twisted tapered cantilever 
beams treated by collocation method", 

Indian Journal of Pure and Applied Physics, 
pp# 459-462, Vol, 10, June 1972, 

1 2* RAO J * S, , 

"Coupled vibxations of turboraachlne blades", 

Shock and Vibration Bulletin, 

pp, 107-125, Bull, 47, Part 2, September 1977, 

13, GUPTA R,S, and RAO S.S., 

"Finite element eigenvalue analysis of tapered and 
twisted Timoshenko beams", 

Journal of Sound and Vibration, 
pp. 187-200, Vol. 56, 1978, 

14, BAJAJ G,R, and tJHOOPAR B.t., ^ 

"Free vibrations of packetted blades — coupled bending 
torsion modes", 

pp# 123-135, 

15, KUtKARNl S,V,, RAO J.S, and SUBRAHMANYAM K.B.- 
"Coupled bending-bending vibrations of pre-twisted 
cantilever blading allomng for shear deflection and 
rotary inertia by the Relssner method". 

International Journal of Mechanical Science, 

Pergamon Press, pp* 517-530, Vol, 23, No, 9, 1981, 

16* KAZA K.R.V. and SUBRAHMANYAM K, B. , ^ ^ 

"Vibration analysis of rotating turbomachlnery blades 
by an Improved finite difference method", . ^ ^ 

International Journal for Numerical Methods in Engineering, 
pp. 1871-1886, Vol* 21, 1985, 

17 , KIEiB H*E,, LEISSA A.l, and MACBAIN J.C,, 

"Vibrations of twisted cantilever plates — A comparison 
of theoretical results", , . r. , « 

International Jou3Ciial for Numerical Methods in Engineering, 
pp. 1365-1380, Vol* 21, ml 8, August 



65 


18, 


19, 


20 , 


21 , 


SAURABH KUMAR, 

"Coupled vibrations of turbine blades by finite element 
method", 

M.Tech, Thesis, IIT KANPUR, 1982* 

TIfOSHENKO S,P, and GOODIER J.N., 

"Theory of Elasticity", 

McGraw Hill Book Co, - Singapore, Third edition, 1982* 
TIKOSHEMCO S*P,, 

"Strength of materials. Part II s Advanced theory and 
problems", 

D, Van Nostrand Co,, Inc*, Third edition, 1956, 

TIMOSHENKO S.P* and WIONOWSKY-KRIEGER S., 

"Theory of Plates and Shells", 

McQraw Hill Book Co,, Singapore, Second edition, 1970. 


22. BROCK A*E. and HOUSTON E.U, 

"Aerodynamic* for Engineering Students", 

Edward Arnold (Publishers) Ltd., 1970* 

23, CHURCH EDWIN F, JR*, 

McQraw Hill Book* Co,, Inc#, Third edition, 1950* 


24, LEE JOHN F*, 

"Theory and Design of Steam and Gas Turbines", 
McQraw Hill Book Co,, Inc,, 1954, 

25. SCHWARZ H.R., RUTISHAUSER H* and STIEFEL E., 
"Numerical Analysis of Symmetric Matrices", 
Prentice-Hall, Inc., 1973. 


26. SOEDEL W,, 

"Vibrations of Shells and Plates", 
Marcel Dicker Inc*, 19© 1, 


27. BREBBIA C.A. and CONNOR J,J., 

"Fundamentals of Finite Elenaent Techniques", 
John Wiley and Sons, New Y^fk-Toronto, 1974, 


28. 


COOK R.D. . ^ . . , . » 
"Concepts and Applications of Finite Elemait Analysis", 
John Wiley and Sons, Inc*, Second edition, 1981, 


29. ZIENCIWia 0,C. _ 

"The Finite Element Method". 

Tata McQraw Hill Book Co, Ltd*, 1979, 


30, NAG MARK 10 User's Manual* 



66 


APPENDIX I 


THE ZIKOVSKY TRANSFORMATION 22 


^ t + 3- 


The Zukovsky transformation 22 

2 

^ where t « x + iy and S’ * % 4- (I*t) 


is a conformal mapping that transforms a circle in the Z'-plane 
into an aeirofoil section in the ^-*plane as shown in the 
Figure 2.3. 

For any point p on the circle in the z plane, 
r » q e cos© + s sin© + t cos^ (I #2) 

For small values of e and , 


(1*3) 


(£ + M) 2 and 

M ^ 

(| ^ |) 2©(1 + cos©) + 2psin© 

Substituting z » r e^® in the Zukovsky transformation# one 
gets 

m 2q cos© and 

2 

y » 2q e( 1 cos©)sin© + 2q sin © 

Using (I*4)f thickness# the area# the area uraments 
of inertia and the product of Inertia can be obtained att 

h •» 4q e sin©(l + cos©) 



A 


4ti e 




'-Tl 


q eCSe'" 4 5 


r,Ti 




. 4 

4-n: q e 


2ti q e P 


(1.5) 


when p is very small compared to e, the inclination of the 
principal axes 1 and 2 v.dth respect to the ^ - r, axes is 
quite small, usually of the order of 2®. Then 


^11 

ro 

Lim 

p -♦ 0 

tl 

4KA 

A 4 3 

6n q e 

^22 

rsJ 

Lim 

P - 0 

^ViT) - 

. 4 

4ti q e 

For 

this 

case 

, the 

expression for q and e are 


( 1 . 6 ) 


a 


(i 


(21^^/3122) 


1/2 


( 1 . 7 ) 


The chord length of the blade comes out to be 
b = 4q : 


4(l22/A)1/2 


( 1 . 8 ) 


The maximum thickness is obtsined by maximising the 
expression h^ of equation (l.5): 

h = 3V3 q e = 3 ( 21 , U.9) 
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appendix II 

householder* s decomposition 


Normally a non*-slngiilar symmetric matrix P is 
decomposed as 

LpJ * luj^ Lq] [u2 (II. 1) 

where [n] is orthogonal and [q] is positive definite. However, 
instead of expressing this transformation as a product invol- 
ving two matrices, one can express it in terms of a matrix 
and a column vector by expressing [uJ as 

L'fl = [l]-2|w][wU <11. 2) 

where W is a normalised coltimn vector. 

It can be found out by the relation 


» 0 , 

W,- U(i 

jCh diC 9 


w,. 



for k«s3, 4 **. n 


vdiere 


* * ^ti jf ^ ^ "*^ij^ (lX.3) 

and the sign in the expression for Wg and Wj^ can be chosen to 
be that of 
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Using (ll*i) ayni th© decomposition for P 

can be expressed is 

Lp] « [q] . 2|W} [ wF [_q 3 ^ 2[Q] {w] {wf 4 

4{wl; 

Advantage® f 

i ) For the higher order matrix P , less computer memory 
locations are requlye<j[ as in lieu of a matrix U of 

the order n x n|t a poitan vector W of n components 
is utilized* 

ii ) Computationti economy and numerical accuracy can be 
achieved fifflultaneouaiy. 
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